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ABSTRACT Recent studies have shown that individuals with colorectal cancer have an altered gut microbiome compared to 
healthy controls. It remains unclear whether these differences are a response to tumorigenesis or actively drive tumorigenesis. 
To determine the role of the gut microbiome in the development of colorectal cancer, we characterized the gut microbiome in a 
murine model of inflammation-associated colorectal cancer that mirrors what is seen in humans. We followed the development 
of an abnormal microbial community structure associated with inflammation and tumorigenesis in the colon. Tumor-bearing 
mice showed enrichment in operational taxonomic units (OTUs) affiliated with members of the Bacteroides, Odoribacter, and 
Akkermansia genera and decreases in OTUs affiliated with members of the Prevotellaceae and Porphyromonadaceae families. 
Conventionalization of germfree mice with microbiota from tumor-bearing mice significantly increased tumorigenesis in the 
colon compared to that for animals colonized with a healthy gut microbiome from untreated mice. Furthermore, at the end of 
the model, germfree mice colonized with microbiota from tumor-bearing mice harbored a higher relative abundance of popula- 
tions associated with tumor formation in conventional animals. Manipulation of the gut microbiome with antibiotics resulted in 
a dramatic decrease in both the number and size of tumors. Our results demonstrate that changes in the gut microbiome associ- 
ated with inflammation and tumorigenesis directly contribute to tumorigenesis and suggest that interventions affecting the 
composition of the microbiome may be a strategy to prevent the development of colon cancer. 

IMPORTANCE The trillions of bacteria that live in the gut, known collectively as the gut microbiome, are important for normal 
functioning of the intestine. There is now growing evidence that disruptive changes in the gut microbiome are strongly associ- 
ated with the development colorectal cancer. However, how the gut microbiome changes with time during tumorigenesis and 
whether these changes directly contribute to disease have not been determined. We demonstrate using a mouse model of 
inflammation-driven colon cancer that there are dramatic, continual alterations in the microbiome during the development of 
tumors, which are directly responsible for tumor development. Our results suggest that interventions that target these changes 
in the microbiome may be an effective strategy for preventing the development of colorectal cancer. 
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Colorectal cancer (CRC) is one of the most commonly diag- 
nosed malignancies worldwide, resulting in over a half- 
million deaths annually (1). Significant risk factors for CRC in- 
clude diets rich in red and processed meat, alcohol consumption, 
and chronic inflammation of the gastrointestinal tract (2-5). Each 
of these factors is closely associated with changes in composition 
and function of the complex community of microorganisms that 
inhabits our gastrointestinal tract. This community, known as the 
gut microbiome, promotes various physiological functions that 
are associated with cancer, including cell proliferation, angiogen- 
esis, and apoptosis (6-9). Therefore, we hypothesized that the 
composition, structure, and functional capacity of the gut micro- 
biome all directly affect tumor development in the colon. 

Several recent studies have addressed this hypothesis by char- 
acterizing the composition of the gut microbiome associated with 



patients with CRC (10-16). Using culture-independent ap- 
proaches, each of these studies observed a significant shift in the 
composition of the gut microbiome in patients with CRC com- 
pared to that in healthy controls. This phenomenon, referred to as 
dysbiosis, can be observed in both the luminal microbiome from 
feces and the mucosa-associated microbiome from tumor biopsy 
specimens. Interestingly, each of these studies obtained conflict- 
ing results regarding the composition and structure of the CRC- 
associated microbial community. Furthermore, there are no bac- 
terial populations that have consistently been identified across 
each study that can be attributed to the development or presence 
of CRC. These data clearly show an association between abnor- 
malities in the gut microbiome and CRC; however, the conflicting 
results point out the need for a mechanistic understanding of the 
role of the gut microbiome in this process. 
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FIG 1 Inflammation-induced tumorigenesis is commensal dependent. (A) Mice were injected with azoxymethane (AOM) on day 1, followed by 3 subsequent 
rounds of water-administered 2% DSS. Colons were harvested 73 days after AOM, and tumors were grossly counted. Black wedges indicate fecal samples used 
for gut microbiome analysis (n = 12). (B) Representative mice were euthanized following each round of DSS to identify macroscopic tumors (n = 5 for each time 
point). An antibiotic cocktail of metronidazole, streptomycin, and vancomycin was administered in the drinking water of a separate cohort of mice for the 
duration of the model (n = 9). Statistical analysis was performed using a two-tailed Student's t test. *, P < 0.01. (C) Representative images of tumors in the distal 
colon of conventional mice treated with AOM/DSS (« = 12) and mice treated with an antibiotic cocktail and AOM/DSS (n = 9). Error bars represent ±SEM. 



The combination of factors that could lead to dysbiosis is com- 
plex and not well understood. In addition, the effect of the devel- 
opment of this abnormal community on colon tumorigenesis re- 
mains unclear. Recent evidence suggests that certain strains of 
Bacteroides fragilis and Escherichia coli can directly affect tumor 
development in the colon through the production of virulence 
factors (e.g., toxins and gene products) (17, 18). Furthermore, 
bacterial populations that produce the short-chain fatty acid bu- 
tyrate have antitumor effects in the colon by promoting apoptosis 
of colonic cancer cells ( 1 9, 20) . We reason that dysbiosis of the gut 
microbiome leads to both enrichment of cancer-promoting bac- 
terial populations and loss of protective populations. Thus, un- 
derstanding the dynamics changes in the gut microbiome on a 
community-wide scale will be essential for understanding colon 
tumor development. 

The gut microbiome is also likely to contribute to CRC 
through the initiation of inflammation. The link between inflam- 
mation and cancer is well established, and patients with inflam- 
matory bowel diseases, such as ulcerative colitis, are at a greater 
risk of developing CRC in their lifetime. In the case of ulcerative 
colitis, the risk for cancer is related to both the duration and se- 
verity of inflammation, with an increasing rate of 0.5 to 1% per 
year after the first decade (2, 21, 22). Chronic inflammation of the 
colon leads to the production of various inflammatory cytokines 
and reactive oxygen species that work in concert to generate a 
tumor microenvironment that promotes carcinogenesis (21, 23, 



24). It has been suggested that this process is microbe driven, but 
it is unclear how the normally beneficial gut microbiome becomes 
inflammatory. 

To determine the role of the gut microbiome in inflammation 
and colon tumorigenesis, we used a well-established model of 
colitis-associated CRC that recapitulates the progression from 
chronic inflammation to dysplasia and adenocarcinoma in hu- 
mans (25). We characterized the dynamics of the gut microbiome 
in this model and demonstrated that community-wide changes 
promote tumorigenesis in the colon. Our data support a model in 
which epithelial cell mutation and inflammatory perturbations of 
the gut microbiome lead to the development of an abnormal mi- 
crobial community with enhanced tumor-promoting activity. 

RESULTS 

Inflammation-associated colon tumorigenesis. We were able to 
replicate an inflammation-based murine model of tumorigenesis 
in specific-pathogen-free (SPF) C57BL/6 mice (n = 12) using an 
intraperitoneal injection of the chemical carcinogen azoxymeth- 
ane (AOM) followed by three subsequent rounds of water- 
administered 2% dextran sodium sulfate (DSS) treatment (26, 27) 
(Fig. 1A). The mice showed a consistent pattern of weight loss 
following each round of DSS treatment, with the most pro- 
nounced change occurring after the first round of DSS (see 
Fig. S1A in the supplemental material). We did not observe mac- 



2 mBio' mbio.asm.org 



November/December 2013 Volume 4 Issue 6 e00692-13 



Gut Microbiome Modulates Colon Tumorigenesis 



roscopic tumors following the first round of DSS administration; 
however, we did observe increased infiltration by immune cells, 
lytter effsignificant epithelial damage, and submucosal edema (see 
Fig. SIB). In addition, we observed a significant increase in the 
proinflammatory mediators macrophage inflammatory protein 2 
(MIP-2), gamma interferon (IFN-y), tumor necrosis factor alpha 
(TNF-a), interleukin 6 (IL-6), and IL-1/3 (see Fig. S1A). Macro- 
scopic tumors and epithelial hyperplasia were apparent following 
the second round of DSS (Fig. IB; see also Fig. SIB). At the end of 
the model, the cohort had a median of 14.5 tumors per mouse (n 
= 12), the majority of which were greater than 1 mm in diameter 
and located in the distal colon and rectum (see Fig. SIB). These 
results demonstrate that our cohort of AOM/DSS-treated mice 
developed a significant number of colonic tumors with complete 
penetrance that could be detected as early as 7 weeks after AOM 
injection. 

To determine whether tumor incidence and penetrance were 
dependent on the gut microbiome, we treated mice (« = 9) with 
an antibiotic cocktail of metronidazole, vancomycin, and strepto- 
mycin ad libitum for 2 weeks prior to AOM and then throughout 
the model, including the days of AOM injection and throughout 
the DSS treatment and recovery periods. Antibiotic-treated mice 
had significantly fewer tumors in the colon than untreated mice 
(Fisher exact test, P< 0.001) (Fig. IB). Tumors that were present 
in antibiotic-treated mice were also significantly smaller than 
those observed in untreated mice (Student's t test, P = 0.002) 
(Fig. 1C; see also Fig. S5A in the supplemental material). These 
results suggest that specific populations within the microbiome 
were essential for tumorigenesis. To determine whether the rela- 
tive change in bacterial density following antibiotic treatment was 
due to a change in the bacterial load, we performed quantitative 
PCR (qPCR) with the 16S rRNA gene from stool samples of 
antibiotic-treated mice. The number of 16S rRNA gene copies per 
mg of feces was not significantly different from that for untreated 
stool samples (P = 0.21) (see Fig. S2). Combined, these results 
indicate that changes to the structure of the community rather 
than total bacterial numbers affected tumorigenesis. 

Significant shifts in the microbiome are associated with co- 
lon tumorigenesis. To further test the hypothesis that specific 
changes in the microbial community structure were associated 
with inflammation and tumorigenesis, we examined the dynamics 
of the gut microbiome throughout the model using stool samples 
from a subset of the original cohort of conventional mice treated 
with AOM/DSS for Fig. 1 (n = 10). We used the fecal samples 
taken prior to AOM injection as a baseline control for each mouse 
and then took samples following each subsequent round of DSS 
administration (Fig. 1A). Mice showed a significant decrease in 
microbial diversity in the gut microbiome following the first 
round of DSS administration through tumor development (P < 
0.001) (Fig. 2A and B). Ordination of the distances between fecal 
samples showed that at the time of euthanization, tumor-bearing 
mice developed a significantly altered microbiome that clustered 
separately from that in baseline samples taken prior to the first 
round of DSS (Fig. 2C). Further examination of fecal samples 
collected at various time points during the AOM/DSS tumor in- 
duction protocol revealed that significant alterations in the micro- 
biome could be observed as early as the first round of DSS admin- 
istration in 7 of the 10 mice. Each round of DSS treatment resulted 
in a significant change in the structure of the microbiome 
(Fig. 2D). Fecal samples taken from tumor-bearing mice after the 



third round of DSS until the time of euthanization also clustered 
separately from earlier samples. The distances between clusters 
were significantly higher than the distances within clusters 
(Fig. 2D). These clusters were observed using operation taxo- 
nomic unit (OTU) and phylogenetics-based metrics of j3-diversity 
(i.e., 0 Y c an d unweighted or weighted UniFrac) and could be dis- 
tinguished from one another using the Random Forest machine 
learning algorithm (accuracy for each group: baseline, 100%; DSS 
round 1, 72.4%; DSS round 2, 71.9%; DSS round 3, 80.6%). These 
results highlighted the association between a dramatically altered 
microbiome structure and the presence of tumors. 

To determine the effect of inflammation on the microbial 
community independent of tumorigenesis, we treated mice with 
three rounds of DSS without the AOM injection (n = 5). There 
was an initial community shift following the first round of DSS, 
but the subsequent stepwise shifts that occurred in AOM/DSS- 
treated mice were not observed in mice treated with DSS only (see 
Fig. S3 in the supplemental material). Furthermore, we did not 
observe the sustained drop in microbial diversity that was ob- 
served in AOM/DSS-treated animals (see Fig. S3A). These results 
suggest that inflammation alone is not sufficient to cause micro- 
bial community changes. Rather, the synergistic effects of the 
AOM/DSS model are necessary for the development of the altered 
microbiome structure and tumorigenesis. 

We next identified which OTUs were responsible for the dra- 
matic shifts in the microbial community structure during inflam- 
mation and tumorigenesis (Fig. 3). Consistent with our commu- 
nitywide j3-diversity analyses, we observed changes in 37 bacterial 
populations (after excluding OTUs representing <0.5% of the 
community) during the time course of the model relative to those 
in baseline samples prior to treatment. Fecal samples taken after 
the first round of DSS were enriched in the relative abundance of 
OTUs affiliated with members of the genus Bacteroides (OTUs 1 
and 13). We also observed a significant decrease in the relative 
abundances of OTUs associated with members of the genus Pre- 
votella and unclassified genera within the family Porphyromon- 
adaceae. Following the second round of DSS, we observed a fur- 
ther loss of the same Prevotella (OTUs 4 and 5) and 
Porphyromonadaceae (OTUs 7, 12, 15, 22, 31, and 48) and the 
continued enrichment of Bacteroides (OTUs 1 and 13). Samples 
taken from mice following the third round of DSS showed signif- 
icant differences compared to those taken following the first 
round of DSS and from healthy baseline mice (Fig. 3 ) (all P values 
were <0.001 as determined by analysis of molecular variance 
[AMOVA]). Tumor-bearing mice showed enrichment in OTUs 
affiliated with Bacteroides (OTU 1), Odoribacter (OTU 3), and 
Turicibacter (OTU 20) . Additionally, we detected a marked bloom 
of a member of the Erysipelotrichaceae family (OTU 26), which 
was undetectable in all of the mice prior to the second round of 
DSS, when tumors are not evident. Simultaneous with the bloom- 
ing of several bacterial populations, there was a significant de- 
crease in the relative abundance of OTUs associated with mem- 
bers of the genus Prevotella (OTUs 4 and 5) and the family 
Porphyromonadaceae (OTUs 7, 12, 15, 22, 31, and 48). An OTU 
associated with the Bacteroides genus (OTU 13), which bloomed 
during the onset of inflammation, decreased significantly follow- 
ing the third round of DSS. These results strongly suggest that 
both inflammation and tumorigenesis promote gut microbiome 
dysbiosis, as highlighted by major shifts in bacterial populations 
from a wide range of taxonomic groups. 
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FIG 2 Development of a dysbiotic gut microbiome during colon tumorigenesis. Microbiome analysis was performed with fecal samples from 10 representative 
mice; color coding is as indicated in Fig. 1 A. (A) Inverse Simpson's diversity index. (B) Observed community richness estimate. Statistical analysis was performed 
using repeated-measures paired group analysis of variance. (C) Nonmetric multidimensional scaling (NMDS) ordination based on 8 yc distances for all 10 mice 
during the AOM/DSS model. (D) Average 0 YC distance within (black) and between (gray) phases of the model. Error bars represent ±SEM. 



We hypothesized that the variability in tumor burden among 
AOM/DSS-treated mice was associated with variability in the gut 
microbiome between mice (coefficient of variation for tumor bur- 
den = 37.9) (Fig. 1C). We identified an OTU related to an unclas- 
sified genus within the family Porphyromonadaceae (OTU 12) that 
was negatively correlated with tumor burden (Spearman correla- 
tion = —0.73; P < 0.05). The relative abundance of this bacterial 
population decreased with each round of DSS, and this drop in 
abundance was more pronounced in mice with higher tumor bur- 
dens. These results suggest that alterations in the relative abun- 
dances of specific bacterial populations were associated not only 
with the incidence of tumors but also with their prevalence. 

Tumor-associated alterations in the microbiome increase 
tumorigenesis in germfree mice. To determine whether the 
community- wide microbiome changes directly contributed to tu- 
mor incidence in the colon, we conventionalized germfree mice 
with either the healthy microbiome of untreated mice or the mi- 
crobiome of tumor-bearing mice analyzed in Fig. 1 . To ensure that 
mice were repeatedly inoculated and stably colonized, we trans- 
ferred fresh feces and bedding to two groups of germfree mice (« 
= 10/group). One group was housed with the bedding from 
healthy, untreated SPF mice, and a second group was housed with 
bedding from tumor-bearing AOM/DSS-treated mice. To mini- 



mize litter effects, each group was comprised of two cages of 5 
mice collected from separate litters that were randomly assigned 
to each of the cages. Following conventionalization, mice were 
treated with AOM/DSS under germfree conditions, as described 
above (Fig. 1). All bacterial phyla and 90% (62 of 69) of genus- 
level taxa detected in donor samples were detected within the re- 
cipient germfree mice (see Table S3 in the supplemental material), 
which is higher than has been previously reported (28, 29). Fur- 
thermore, 81% of the sequences we obtained from the donor mice 
belonged to OTUs that were found in the recipient germfree mice. 
Mice conventionalized with the microbiome of tumor-bearing 
mice had a 2-fold increase in tumor burden (P = 0.002) relative to 
that of mice conventionalized with a healthy microbiome 
(Fig. 4A). Additionally, tumors from these mice were significantly 
larger than those observed in recipients of a healthy microbiome 
(P = 0.002) (see Fig. S5B). Similar to our results with SPF mice, 
germfree mice conventionalized with the community of tumor- 
bearing mice had a significantly less diverse gut microbiome (P < 
0.001). Using community-wide j8-diversity analyses, we deter- 
mined that conventionalization with these two treatments of bed- 
ding resulted in two distinct microbial community structures 
(AMOVA. P < 0.001) (Fig. 4C). Germfree mice conventionalized 
with the microbiome of tumor-bearing mice showed significant 



4 mBio' mbio.asm.org 



November/December 2013 Volume 4 Issue 6 e00692-13 



Gut Microbiome Modulates Colon Tumorigenesis 



Time 





Change in 
relative abundance 

i , , , n 

-0.05 0 0.05 0.1 



OTU 


01 


(Bacteroides) 


OTU 


02 


(Odoribacter) 


OTU 


13 


(Bacteroides) 


OTU 


20 


(Turicibacter) 


OTU 


35 


(Clostridium) 


OTU 


26 


(Erysipelotrichaceae) 


OTU 


23 


(Alistipes) 


OTU 


40 


(Clostridiaceae) 


OTU 


19 


(Parabacteroides) 


OTU 


42 


(Lachnospiraceae) 


OTU 


24 


(Parabacteroides) 


OTU 


31 


(Porphyromonadaceae) 


OTU 


27 


(Anaeroplasma) 


OTU 


15 


(Porphyromonadaceae) 


OTU 


22 


(Porphyromonadaceae) 


OTU 


49 


(Bacteria) 


OTU 


51 


(Lactobacillus) 


OTU 


09 


(Bacteroidetes) 


OTU 


05 


(Prevotella) 


OTU 


12 


(Porphyromonadaceae) 


OTU 


04 


(Prevotellaceae) 


OTU 


48 


(Porphyromonadaceae) 


OTU 


07 


(Porphyromonadaceae) 



FIG 3 Heat map of OTUs with relative abundances that are significantly different from their relative abundances at the time of AOM administration. The 
average OTU abundance between mice for each OTU was calculated for each time point. The timeline is colored for the following groups: baseline samples (prior 
to AOM), black; following the first round of DSS, blue; following the second round of DSS, green; following the third round of DSS, red. The OTU number and 
taxonomic group based on RDP classification are represented for each row. Repeated-measures paired group analysis of variance was used to identify significantly 
altered OTUs. 



enrichment in the relative abundance of OTUs affiliated with the 
genera Bacteroides (OTU 1) and the family Erysipelotrichaceae 
(OTU 26). Additionally, these germfree mice had significantly 
fewer members of the Porphyromonadaceae (OTU 12) compared 
to germfree mice conventionalized with bedding from healthy 
mice. Finally, germfree mice conventionalized with a healthy mi- 
crobiome successfully recapitulated the community dynamics 
seen in conventional mice during tumorigenesis. We observed 
significant changes in 34 OTUs following the AOM/DSS model. 
Similar to tumor-bearing conventional mice, germfree tumor- 
bearing mice showed enrichment in OTUs affiliated with mem- 
bers of the Bacteroides (OTU 1), Odoribacter (OTU 3), Turicibac- 
ter (OTU 20), and a bloom in the Erysipelotrichaceae (OTU 26), 
which was undetectable before AOM administration. There was 
also a significant decrease in the relative abundance of OTUs as- 
sociated with members of the genus Prevotella (OTUs 4 and 5) and 
the family Porphyromonadaceae (OTUs 6, 7, and 12). These results 
demonstrate that alterations to the gut microbiome that were as- 
sociated with chronic inflammation and tumorigenesis in SPF 
mice were transmitted to germfree mice and can exacerbate colon 
tumorigenesis. 



DISCUSSION 

In the present study, we established a causal role for the gut mi- 
crobiome in exacerbating tumor formation in an inflammation- 
based model of tumorigenesis. Manipulation of the microbiome 
using antibiotics reduced tumor formation, which highlighted the 
importance of bacterium-driven factors in tumorigenesis. We 
demonstrated dynamic changes in the microbial community 
structure associated with dysbiosis, which occurs prior to the first 
signs of macroscopic tumor formation. We established the syner- 
gistic effect of AOM- and DSS-induced inflammation and tumor- 
igenesis in driving microbial community changes that occur in a 
stepwise fashion. Finally, transfer of microbiota from tumor- 
bearing mice into germfree mice significantly increased the num- 
ber and size of tumors compared to those in germfree mice inoc- 
ulated with healthy microbiota. Our experiments also 
demonstrated dramatic shifts in the relative abundances of bacte- 
rial populations, including those related to the genus Bacteroides, 
that were associated with increased tumorigenesis. 

Several recent studies have compared the gut microbiome of 
patients with CRC to that of healthy controls (10-16). These stud- 
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FIG 4 Tumor-associated gut microbiome alterations exacerbate tumorigenesis in germfree mice. (A) Number of tumors observed at the end of the model when 
germfree mice were colonized using bedding from healthy mice (Healthy community) or mice with tumors (Dysbiotic community). (B) Representative images 
of tumors in the distal colon of mice conventionalized with a healthy microbiome (n = 10) or the microbiome of tumor-bearing mice (n = 9). (C) NMDS 
ordination based on 0 yc distances for all 19 mice following conventionalization with a healthy microbiome (Healthy community) or the microbiome of 
tumor-bearing mice (Dysbiotic community). Error bars represent ±SEM. 



ies have consistently demonstrated significant differences in the 
microbial community structure of patients with CRC, but each 
study has disagreed in terms of the specific gut microbiome com- 
position and profile associated with CRC. The inability to identify 
a consensus community profile or etiological agent is likely due to 
the large variation in the structure of the microbiome across indi- 
viduals and the improbability of there being one community pro- 
file or bacterial population that is associated with all CRCs. We 
were able to reduce the interindividual variation and the diversity 
of cancer types using a murine model of inflammation-induced 
CRC. Unlike the human cross-sectional studies, we demonstrated 
dynamic changes in the microbiome during the development of 
inflammation and tumorigenesis and that these changes directly 
cause disease. 

Based on this study, the gut microbiome complements the ac- 
tivity of AOM and DSS to cause tumorigenesis, but the underlying 
mechanisms driving microbially mediated tumorigenesis ob- 
served remain to be elucidated. Although a number of bacterial 
populations have altered relative abundances throughout the 
model, it is as yet unclear whether there is an increase in bacterial 
populations that induce inflammation or a loss of populations 
that produce anti-inflammatory signals and help maintain im- 
mune homeostasis in the gut. Regardless, an increasingly inflam- 
matory environment would generate a self-reinforcing patho- 



genic cascade between the gut microbiome and the host, fostering 
the development of cancer through the development of, for exam- 
ple, genotoxic reactive oxygen species and protumor inflamma- 
tory mediators (e.g., TNF-a, IL-6, IL-1/3, and IL-23). In addition 
to the role of the gut microbiome in inflammation, changes me- 
diated by chronic inflammation and tumorigenesis could lead to 
the enrichment of bacterial populations (30) that have a direct role 
in tumor development through the production of metabolites, 
antigens, virulence factors, and other potential tumor-promoting 
gene products. A recent study by Arthur et al. (17) demonstrated 
that colonic inflammation in the IL-10-deficient mouse impacts 
the composition of the gut microbiome, leading to an enrichment 
of tumor-promoting E. coli strains. Although we did not detect 
any significant changes in populations related to the genus Esche- 
richia, it is likely that the microbial community alterations we 
observed in our tumor model are enriched with populations that 
fill a similar role. Specifically, marked increases in Bacteroides spe- 
cies in our study may contribute to tumorigenesis. Human com- 
mensals belonging to the genus Bacteroides, specifically entero- 
toxigenic B. fragilis (ETBF), have been associated with 
inflammation and CRC (18,31). ETBF has been shown to strongly 
induce colonic tumors in multiple-intestinal-neoplasia mice 
through secretion of a metalloprotease toxin, and certain strains 
are thought to contribute to CRC risk in humans. We did not 
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detect ETBF in the murine gut microbiome (data not shown), but 
it is possible that similar processes are occurring during tumori- 
genesis in mice. 

Chronic inflammation and tumorigenesis are also likely to lead 
to the loss of members of the gut microbiome that are important 
for maintaining epithelial health and immune homeostasis. In this 
study, we observed a dramatic decrease in OTUs from unclassified 
genera within the family Porphyromonadaceae. We hypothesize 
that these bacterial populations serve a protective role and are 
important mediators of gut health in the murine gut microbiome. 
One mechanism of protection could be through the fermentation 
of complex carbohydrates (e.g., fiber) into short-chain fatty acids 
(SCFA), such as butyrate. Butyrate reduces inflammation (32) and 
inhibits growth and induces apoptosis in cancer cells (19, 20). 
Therefore, loss of butyrate-producing populations in the gut 
could increase both inflammation and tumorigenesis. This is sup- 
ported by extensive epidemiological data that demonstrate a link 
between diets high in fiber and a decreased CRC risk (33). Fur- 
thermore, recent studies have shown that individuals who con- 
sume low-fiber diets or are diagnosed with CRC have a lower level 
of SCFAs in their feces (34). It is also possible that members of the 
family Porphyromonadaceae are important mediators of anti- 
inflammatory signals in the gut. A loss of such anti-inflammatory 
populations would lead to a dramatic intensification of inflamma- 
tion in the gut during DSS-induced colitis and a marked increase 
in tumor-promoting signals. 

It is important to note that the gut microbiome is an extremely 
complex and diverse community, and therefore it is unlikely that 
a single bacterial population is responsible for driving tumorigen- 
esis or that one CRC-associated microbiome can be found in all 
CRC patients. Rather, as our data suggest, a community-wide ef- 
fect involving the gain and loss of bacterial populations and gen- 
eral metabolic functions likely plays a critical role in CRC devel- 
opment. As we demonstrated in this study, changes in the entire 
gut microbiome can dramatically alter the tumor burden, and 
identifying the mechanisms behind this phenomenon will be crit- 
ical for addressing how the microbiome can be altered therapeu- 
tically to reduce colon tumorigenesis. 

MATERIALS and METHODS 

Animals and animal care. Studies were conducted on adult (8 to 12 weeks 
old) age-matched male C57BL/6 mice that were bred and maintained 
under SPF or germfree conditions as specified above. Both SPF and germ- 
free mice were fed the same autoclaved chow diet. All animal experiments 
were approved by the University Committee on Use and Care of Animals 
at the University of Michigan. 

Inflammation-induced colon tumorigenesis. Eight- to twelve-week- 
old mice received a single intraperitoneal (i.p.) injection of azoxymethane 
(10 mg/kg of body weight). Water containing 2% DSS was administered 
to mice beginning on day 5 for 5 days, followed by 16 days of water. This 
was repeated twice for a total of 3 rounds of DSS (26, 27). Mice were 
euthanized on days 14, 24, 38, and 45 for intermediate time point analysis. 
The remaining mice were euthanized 3 weeks after the third round of DSS 
administration for tumor counting. 

Histological analysis. At necropsy, all colons were harvested, flushed 
of luminal contents, and cut open longitudinally to count and measure 
tumors. The largest dimension of each tumor was measured with calipers. 
Tumors were categorized based on size (<1 mm, 1 to 2 mm, or >2 mm). 
Colons were then "jelly rolled," fixed in Carnoy's solution, and embedded 
in paraffin. Five-micrometer sections were used for hematoxylin-and- 
eosin (H&E) staining, and slides were analyzed under magnification 
X100. 



RNA isolation and cytokine analysis. Distal colon tissue was homog- 
enized and total RNA was isolated using the Nucleospin RNA kit 
(Macherey-Nagel). cDNA was synthesized using an iScript kit (Bio-Rad). 
and the cDNA was then used for quantitative PCR (qPCR) using a SYBR 
green expression assay (Applied Biosystems). 

Antibiotic treatment. Mice were treated with an antibiotic cocktail of 
metronidazole (0.75 g/liter), vancomycin (0.5 g/liter), and streptomycin 
(2 g/liter) in their drinking water for 2 weeks prior to and throughout the 
duration of AOM/DSS administration. 

Germfree conventionalization. Eight-weeks-old C57BL/6 germfree 
male mice were used. Fresh feces and bedding were collected from un- 
treated and AOM/DSS-treated tumor-bearing mice and immediately 
transferred to cages of germfree mice 2 weeks prior to AOM injection to 
allow stable colonization. Germfree mice were divided into two treatment 
groups, one group receiving bedding from untreated, healthy mice and 
the other group receiving bedding from AOM/DSS-treated tumor- 
bearing mice. To ensure that there were no cage effects, each treatment 
group was comprised of two cages of mice. The mice were obtained from 
separate litters and randomly assigned to the four cages. Mice were then 
treated with AOM/DSS to induce tumors as described above. Three weeks 
after the last round of DSS, mice were euthanized and colons were har- 
vested as described above. 

DNA extraction. Fecal samples were collected daily from the mice 
throughout the AOM/DSS protocol and immediately frozen for storage at 
— 20° C. We selected 12 fecal samples distributed over the 73-day timeline 
of the AOM/DSS model for 10 representative mice. Microbial genomic 
DNA was extracted using the PowerSoil-htp 96 Well Soil DNA isolation 
kit (Mo Bio laboratories) using an EpMotion 5075 pipetting system. 

16S rRNA gene sequencing and curation. The V35 region of the 16S 
rRNA gene from each sample was amplified and sequenced using the 454 
Titanium sequencing platform at the Baylor College of Medicine Human 
Genome Sequencing Center as described elsewhere (http://www.mothur 
.org/wiki/454_SOP). We curated our sequences as described previously 
using the mothur software package (35-37). Briefly, we denoised se- 
quences using the PyroNoise algorithm after trimming each flowgram to 
450 flows (38), aligned the resulting sequences to a reference alignment 
derived from the SILVA 16S rRNA sequence database (39), and removed 
sequences that were flagged as possible chimeras by UCHIME (50) or that 
did not align to the V35 region. After curating the sequence data, we 
obtained between 6 and 10,742 sequences (median = 5,681), with a me- 
dian length of 253 bp. To minimize biased effects of uneven sampling, we 
rarefied to 1,800 sequences per sample. Seven samples either did not pass 
through sequence curation or had less than 1,800 sequences and were 
therefore not used for further analysis. Parallel sequencing of a mock 
community allowed us to measure a median error rate of 0.06%. 

Analysis of the microbiome. Sequences were clustered into OTUs 
based on a 3% distance cutoff using the average-neighbor algorithm. All 
sequences were classified using the RDP training set, version 9 (http: 
//sourceforge.net/projects/rdp-classifier/), and OTUs were assigned a 
classification based on which taxonomy had the majority consensus of 
sequences within a given OTU using a naive Bayesian classifier (40). Mi- 
crobial diversity was calculated using the inverse Simpson index (41) and 
the observed number of OTUs. To calculate j3-diversity, we used the 0 YC 
distance metric with OTU frequency data (42), and we calculated UniFrac 
statistics (43) using neighbor-joining phylogenetic trees generated using 
the nonheuristic neighbor-joining algorithm implemented in the soft- 
ware program Clearcut (44). Analysis of molecular variance (AMOVA) 
was performed to determine significance between the community struc- 
tures of different groups of samples based on 8 YC and UniFrac distance 
matrices (45). To identify OTUs important for driving differences be- 
tween groups (baseline, after DSS round 1, after DSS round 2, and after 
DSS round 3), we used a repeated-measure paired treatment analysis of 
variance for each OTU and corrected for multiple comparisons using an 
experiment-wise error rate of 0.01 (46). Additionally, we identified fea- 
tures (OTUs) important for each group using the machine learning algo- 
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rithm random forest as implemented in the software environment R 
(http://CRAN.R-project.org) (47). The abundance-based Jaccard dissim- 
ilarity index indicates the fraction of all sequences that affiliate with OTUs 
that are shared between two communities and was used to calculate the 
fraction of OTUs that were shared between donor samples and germfree 
recipient samples (48). All sff files and the MIMARKS spreadsheet are 
available at http://www.mothur.org/aomdss_dynamics/. 

16S rRNA qPCR analysis. Relative bacterial loads in stool samples 
were quantified by qPCR analysis of bacterial genomic DNA using Kapa 
SYBR-fast master mix (Kapa biosciences) and universal 16S rRNA gene 
primers (F, ACTCCTACGGGAGGCAGCAGT; R, ATTACCGCGGCTG 
CTGGC) (49). Samples were normalized to fecal mass, and the relative 
fold change was determined using untreated stool samples for each repli- 
cate mouse (n = 5). Note that qPCR measures the relative fold change of 
the 16S gene copy number, not actual bacterial numbers. 

SUPPLEMENTAL MATERIAL 

Supplemental material for this article may be found at http://mbio.asm.org 
/lookup/suppl/doi: 10.11 28/mBio.00692- 1 3/-/DCSupplemental. 

Figure SI, PDF file, 0.5 MB. 

Figure S2, PDF file, 0.1 MB. 

Figure S3, PDF file, 0.1 MB. 

Figure S4, PDF file, 0.1 MB. 

Figure S5, PDF file, 0.1 MB. 

Table SI, DOCX file, 0.1 MB. 

Table S2, DOCX file, 0.1 MB. 

Table S3, DOCX file, 0.2 MB. 
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